Impurity effects in quantum dots: Towards quantitative modeling 
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We have studied the single-electron transport spectrum of a quantum dot in GaAs/AlGaAs res- 
onant tunneling device. The measured spectrum has irregularities indicating a broken circular 
symmetry. We model the system with an external potential consisting of a parabolic confinement 
and a negatively charged Coulombic impurity placed in the vicinity of the quantum dot. The model 
leads to a good agreement between the calculated single-electron eigenenergies and the experimental 
spectrum. Furthermore, we use the spin-density-functional theory to study the energies and angular 
momenta when the system contains many interacting electrons. In the high magnetic field regime 
the increasing electron number is shown to reduce the distortion induced by the impurity. 



I. INTRODUCTION 

The tunability in size, shape, and electron number of 
semiconductor quantum dots (QD) provides numerous 
technological applications as well as interesting many- 
electron physicsi In actual QD devices, the effects in- 
duced by impurities or donor scattering centers may 
be remarkable. In most cases, irregularities in samples 
have only an indirect influence on the many-body struc- 
ture, complicating the identification of the origin behind 
the peculiar behavior in the measured characteristics of 
QD's. 

A clean quantum dot typically shows single-electron 
energy levels reminiscent of the well-known Fock-Darwin 
(FD) energy spectrum corresponding to a parabolic con- 
fining potential,- Adding external impurities into the QD 
breaks the circular symmetry of the system, leading to 
avoided crossings and liftings of the degeneracies in the 
single-electron energy spectrum. This was demonstrated 
by Halonen et al.^ who studied theoretically QD's dis- 
torted by repulsive Gaussian scattering centers. How- 
ever, even if clear traces of the FD spectrum have been 
obtained experimentally in both lateral^ and vertical^ 
quantum dots, there is, to the best of our knowledge, 
no direct experimental evidence of repulsive impurities 
present in QD structures. Instead, states bound to 
hydrogenic impurities, probably arising from Si dopant 
atoms in the GaAs quantum well, were found already 
by Ashoori and co-workersi in their pioneering single- 
electron tunneling experiment. These impurities have 
been suggested to be sources of pair-tunneling states, the- 
oretically analyzed with a superimposed attractive l/r- 
type potentiali^ 

Theoretically, the distortion of the circular symmetry 
makes the many-electron problem particularly complex 
to solve, especially in the presence of an external mag- 
netic field. In the above-mentioned study, Halonen et al.^ 
applied exact diagonalization up to three electrons and 
focused on the effects of impurities on the energy lev- 
els and optical absorption spectra. Recently, Giiglii and 
co-workergifi performed diffusion quantum Monte Carlo 
(DMC) calculations on QD's distorted by randomly dis- 



tributed Gaussian scatterers and studied the energetics 
up to ten electrons. They found that in these systems 
the transitions between the many-body states are con- 
siderably less pronounced than in clean dots. Hirose and 
Wingreen^^ have used the spin-density- functional theory 
(SDFT) to examine the energies and spin states in dis- 
ordered QD's as a function of the interaction strength in 
zero magnetic field. Besides additional scatterers, non- 
circular QD's have attracted general interest in connec- 
tion with the chaotic propertiesi^ or the behavior in the 
high magnetic field limit ^ 

In this paper we present a measured single-electron 
transport spectrum where avoided crossings and lifted 
degeneracies are clearly observable. We reproduce the 
spectrum with an appropriate model potential, showing 
that the unexpected effects in the spectrum result from 
a negatively charged Coulombic impurity located near 
the QD. The many-electron properties studied by the 
SDFT reflect the strongly distorted single-electron spec- 
trum. The variation of the impurity location shows the 
stability of the maximum-density droplet (MDD) and the 
screening of the impurity by electrons. 

The outline of this paper is as follows. In Sec. Hll we 
briefly describe the fabrication of the sample and report 
the transport measurement. In Sec. IIIII the theoretical 
model describing the physical system is given and the 
single-electron calculations are compared to the experi- 
ment. In Sec. lIVI the many-electron properties, i.e., chem- 
ical potentials, MDD stability, and total magnetization 
are studied with the SDFT. The paper is summarized in 
Sec. El 



II. EXPERIMENT 

The heterostructure consists of a 10 nm wide GaAs 
quantum well sandwiched between two Alo.3Gao.7As- 
tunneling barriers of 5 and 8 nm, see Fig.^ The contacts 
are formed by 0.5 fim thick GaAs layers highly doped 
with Si up to 4 X 10^^ cm"'^ and separated from the active 
region by 7 nm thin spacer layers of undoped GaAs. Our 
sample was defined as a mesa of 40 ^m size. We carried 
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FIG. 1: Sketch of the heterostructure of our sample. 
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FIG. 2: Top; G(V,B) plot of the transport spectrum of our 
sample. Bottom: I-V characteristics for _B = T. 



out direct-current measurements of the current-voltage 
(I-V) characteristics in a He'^-refrigerator at 350 mK base 
temperature in magnetic fields up to 14 T. 

Figure El shows the resulting transport spectrum of a 
quantum dot formed in a local potential minimum. The 
black lines 'Vn,i correspond to high differential conduc- 
tance G = d//dV^. They trace the position of the single- 
electron energy states E^^i of the spectrum according to 
a relation 



(1) 



where the energy-voltage conversion factor a equals 0.4, 
determined from measurements of the broadening of the 
step edge with temperature, and the onset voltage is fit- 
ted to be Vq = 172 mV. The energy levels Vn^i in the 
transport spectrum can be interpreted as single-particle 
levels of a local, presumably a grown-induced potential 
minimum in the GaAs quantum well of our device. The 



lowest FD band is clearly visible, and one can also ob- 
serve higher excited states. In contrast to ordinary FD 
energy levels, we are able to observe broken energy de- 
generacies at B = T and strong anticrossing effects in 
the spectrum. 



III. MODELING THE QUANTUM DOT 

We expect the quantum well confined in the GaAs layer 
to have a negligible degree of freedom for electrons in the 
vertical direction. Our model system is thus strictly two- 
dimensional and defined to be located on the xy plane. 
The single-electron Hamiltonian is written as 



Mr) 
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[p + eA(r)]' + Konf (r) + V^mp(r), (2) 



where we use the effective-mass approximation with 
TO* = 0.067 TOg, which is the typical value for electrons 
moving in GaAs. In a symmetric gauge the vector poten- 
tial reads as A — y,x,0), giving the external mag- 
netic field B — Bz perpendicular to the QD plane. The 
Zeeman energy is omitted in Eq. Q since the spin split- 
tingl^ is not visible in the energy levels shown in Fig. [3 
for the magnetic fields applied. 

The confining potential Vconf(i") is expected to be 
parabolic near the center of the dot. However, we 
soften the edges of the dot by changing the sign of the 
paraboloid at a certain cusp radius Vc, giving 



K:onf(r) 



r <rc 



(3) 

where the parameter s defines the strength of the round- 
ing term. As shown below, the softening of the confine- 
ment is crucial in obtaining a good agreement with the 
experimental energy spectrum. 

We expect the impurity to be described by a negatively 
charged particle located in the vicinity of the quantum 
well. The impurity potential can thus be written in a 
Coulombic form as 



Vimp(r) 



R) 



(4) 



Vimp , and a sketch of 



where q is the (negative) charge of the impurity particle, 
K describes the "dielectricity" between the impurity and 
the electrons in the QD, and R and d are the lateral and 
vertical distances of the impurity from the QD center, re- 
spectively. Figure|31shows the total external confinement 
of the model system, Voxt = K;onf 
the expected configuration. 

To calculate the single-electron spectrum, we solve the 
discretized eigenvalue problem hipi — eiipi [i = 1 ... 12) 
numerically on a 2D point grid using a Rayleigh quotient 
multigrid methodJ-^ Figure 01 shows the resulting spec- 
trum (dashed lines) compared to the experimental data 
(repeated from Fig. The energies are converted to 



FIG. 3: (a) Profile of tfie external potential used in the sim- 
ulation, (b) Sketch of the expected configuration of the QD- 
impurity system. 
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FIG. 4: Measured transport spectrum (repeated from Fig.|5J 
of a GaAs/AlGaAs QD and the calculated single-electron en- 
ergies (dashed lines) corresponding to the model potential 
shown in Fig.|H^a). 



voltages according to Eq. ^ , and the model parameters 
are adjusted (see the discussion below) until the agree- 
ment between the experiment and the model is as good as 
possible. The simulation places the avoided crossings be- 
tween the energy levels very close to the correct positions. 
There are still considerable deviations in the 5th and 6th 
levels but, for example, the 7th level agrees almost per- 
fectly through the magnetic-field regime presented. The 
differences at high fields between the experimental data 
and the simulation result from the shift of the chemical 
potential of the emitter to higher energies with increasing 
magnetic field. 

In calculating the energy spectrum shown in Fig. 0] we 
optimize the potential parameters corresponding to the 
best possible fit to the experimental data. The confine- 
ment is then defined by Hloq = 13.8 meV, Tc = 15.5 nm, 
and s = —0.2, and the impurity parameters are given by 



FIG. 5: Lowest noninteracting single-electron eigenenergies 
calculated for an impurity-containing QD with (thick lines) 
and without (thin lines) the rounding of the edges. The 
dashed lines show the corresponding eigenenergies (lifted by 
-1-14 meV for clarity) for a clean dot. 



q/n — — 2e , R — 14.5 nm, and d — 2 nm. There is nat- 
urally some uncertainty in the parameters. However, for 
this particular sample we find the following characteristic 
features in the model potential: (i) The K-reduced impu- 
rity charge is so large that it indicates a multiple charged 
particle distorting the QD. (ii) The impurity is located 
very close to the QD plane, probably lying in the 10 nm 
thick GaAs layer (see Fig.^. (iii) The confining poten- 
tial is approximately three times larger than the values 
typically used (3 ... 5 meV) for modeling parabolic QD's. 
This is due to the grown-induced formation of the QD in 
the absence of gates around the sample, (iv) For the 
same reason, the confinement becomes softer toward the 
edges of the dot. Hence, the rounding at r > Tc in Vconf 
is required to compress the highest states in agreement 
with the experimental spectrum. 

To clarify the sensitivity of the single-electron spec- 
trum on the shape of the model potential, we compare in 
Fig-Elthe eigenenergies given by the chosen model (thick 
lines) to those of a model QD without the rounded edges 
(thin lines) and to those of a corresponding clean model 
dot without Vimp but with the rounding term (dashed 
lines). The level repulsion is clearly induced by the 
Coulombic impurity that breaks the circular symmetry 
of the QD and mixes the FD levels of different orbital 
quantum numbers. A statistical analysis of the energy- 
level spacings would enlighten the quantum chaotic prop- 
erties^^ of the system but it is not included in this study. 
In the high-field limit, however, the system becomes in- 
tegrable and the eigenstates condense into Landau levels. 

The rounding term in V^onf has the strongest influence 
on the levels with the highest angular momenta, and the 
cusp at Tc induces also a weak decoupling of the degen- 
eracies at S = T. We remark that the eigenenergies for 
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the clean case are lifted in Fig. |Slby 14 meV for clarity. 
The Coulombic impurity in the vicinity of the QD thus 
has a strong effect on the eigenenergies. This tendency 
is also apparent in the many-electron properties studied 
below. 



IV. MANY-ELECTRON PROPERTIES 

Next we study situations that the quantum dot de- 
scribed by the best fitting parameters above contains up 
to six interacting electrons. Even if the many-electron 
case has not yet been experimentally realized for this 
particular QD, we find it important to predict how the 
increasing electron number changes the effects of the im- 
purity on the ground-state properties. 

The problem is now described by the iV-elcctron 
Hamiltonian 



N 



N 



i<j 



47reoe r 



(5) 



where the single-electron part hi is given by Eq. |(2Jl and 
the Zeeman energy is taken into account with g* — —0.22 
for the effective gyromagnetic ratio in GaAs. This has 
been measured to be a realistic value for a similar sys- 
tem. We have e = 12.7 for the dielectric constant in the 
Coulomb interaction between the electrons. 

We apply the SDFT according to the self-consistent 
Kohn-Sham (KS) scheme to obtain the total energies and 
spin densities of the system. The local spin-density ap- 
proximation used for the exchange-correlation energy is 
based on the magnetic-field-independent formulation by 
Attaccalite et ali^ According to our experience, it is the 
most accurate parametrization for finite two-dimensional 
electron systems in the zero-field limit. We solve the 
discretized KS equations in real space without implicit 
symmetry restrictions, which allows us a total freedom 
in shaping the external potential. The numerical process 
of solving the effective single-electron Schrodinger equa- 
tion is accelerated with the Rayleigh quotient multigrid 
methodii^ 

Our earlier calculations for rectangular— and 
parabolioiS, QD's show that our SDFT scheme pro- 
duces energies in a good accordance with the quantum 
Monte Carlo results. We have also noticed that the 
current-spin-density-functional theory (CSDFT) does 
not represent a qualitative improvement over the SDFT 
in small quantum-dot systems. A detailed comparison 
between the SDFT and CSDFT for a six-electron QD 
can be found in Ref . Il9l 



A. Energies 

Figure ini shows the chemial potentials fJ,{N) = E{N) — 
E{N—1), calculated in clean (dashed lines) and impurity- 
containing (solid lines) QD's up to six electrons. Due to 




FIG. 6: Calculated chemical potentials for a clean (dashed 
lines, lifted by -1-14 meV for clarity) and impurity-containing 
(solid lines) quantum dot up to six electrons. The dotted 
lines denote the borders for the full spin polarization. The 
other spin-state domains are also marked above and below 
the curves corresponding to the clean and impurity-containing 
cases, respectively. 



the spin-degeneracy, the chemical potentials shift gen- 
erally in adjacent pairs with increasing B correspond- 
ing to the well-known even-odd effect of experimental 
current peaks^ In the clean case, however, occasional 
pairing of fi{N + 1) and /i(iV — 1) is clearly observ- 
able, e.g., with = 3 and 5, and with A^ = 4 and 
6 at B ^ IT. This is caused by Hund's rule, i.e., 
near a degenerate point it is energetically favorable to 
have parallel spins between the electrons due to the ex- 
change interaction. This leads to partial spin polarization 
(5 = ^ 1) in the clean dot with A^ = 4 and A^ = 6 
at B < 0.84 T and 2.5T < B < 3.5 T, respectively. In 
the impurity-containing dot those states are missing due 
to the avoided level crossings in the single-electron spec- 
trum (see Fig.lSJ. Generally, the impurity smoothes the 
behavior of the chemical potential, reflecting the flatten- 
ing of the single-electron energies studied above. This is 
in agreement with the results by Giiglii et aZ.JJ^ for ran- 
domly distributed impurities. 

As seen in Fig.|Blthe impurity does not affect consider- 
ably the onset for the full spin polarization (dotted lines). 
As A^ increases, however, this point in the clean QD shifts 
to higher B than in the impurity-containing dot. Simul- 
taneously, the fi spacings decrease relatively more rapidly 
in the latter case than in the clean dot. The reason is the 
fact that the impurity pushes the electrons to the soft- 
conflnement region near the edges, whereas the clean dot 
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FIG. 7; Total angular momentum and electron density dis- 
tributions for a six-electron QD as a function of the impurity 
distance from the dot center. The solid and dashed curves 
correspond to the magnetic fields of 20 and 25 T, respectively. 



represents a more compact state with higher addition en- 
ergies. 



B. Angular momentum 

In the description above, the impurity is kept in a fixed 
position at i? = 14.5 nm laterally from the dot center. 
Next, we consider changes in the total angular momen- 
tum as R is varied (while d is kept constant) in mag- 
netic fields corresponding to the MDD state. In a cir- 
cularly symmetric QD, the MDD has the occupation on 
the single-particle states with angular momentum eigen- 
values I = 0,-1, ...,-N+l, giving = -^N{N-1) for 
the total angular momentum^ In the case of a broken 
circular symmetry, the onset for the MDD can be defined 
from the last cusp in the chemical potential indicating 
the full spin polarization (see Fig. ^ . This is consistent 
with the experimental identification of the MDD phase 
by measuring the magnetic field evolution of the Coulomb 
blockade peaks. Since the angular momentum I is not 
a good quantum number for a non-circular confinement, 
we compute Lz in the following as a sum of the expec- 
tation values for the single-electron angular momentum 
operator, = —ih[x{d/dy) — y{d/dx)]. The summation 
is over all the occupied KS states. 

In Fig. [3 we present how the six-electron MDD at 
S = 20 T is affected by the impurity shifted through 
the whole QD region. Except for the kink at i? ~ 8 nm, 
corresponding to the radius at which the impurity pene- 
trates the edge of the QD, \Lz\ increases smoothly from 
15 to 21, i.e., from the MDD to a quantum ring with 
I = — 1,...,— 6 states occupied. The stability is also vi- 
sualized in the density insets of Fig. corresponding to 
different impurity locations. The dashed line shows the 
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FIG. 8: Magnetization of the (fully-polarized) state of a six- 
electron quantum dot at the impurity distances R = (plus 
markers) and 5 nm (cross markers) with (solid curves) and 
without (dashed curves) the electron-electron interactions. 
The effective Bohr magneton /ig = eh/2m* . 



equivalent evolution at i? = 25 T, still corresponding to 
the MDD region in the clean dot. As i? ^ 0, however, 
we find a transition to the next quantum-ring state with 
I = —2,...,— 7 states occupied. The rounded edges of 
the QD clearly make the ring-like states sensitive to the 
magnetic field. In contrast, the MDD remains as a rather 
compact state spread across the highly-confined central 
region as the magnetic field is varied. 



C. Magnetization 

Finally, we investigate the tendency of the many- 
electron structure to screen the external impurity. For 
that purpose, we consider the total magnetization, de- 
fined at zero temperature as M — —OEtot/dB. In Fig. |S| 
we plot M as a function of B for six-electron, fully- 
polarized QD's in both noninteracting and interacting 
cases and with the impurity locations R = and 5 nm 
laterally from the QD center. The noninteracting en- 
ergy is computed directly as a sum of the six lowest 
single-electron energies. As the impurity is shifted off the 
dot center, the lowest single-electron states decouple (see 
Fig. 1 in Ref. 0), leading to a considerably smaller \M\ 
and disappearance of the transition at i? 23 T in the 
noninteracting systems. On the other hand, the interact- 
ing systems show frequent oscillations due to the entan- 
glement of the energy levels, i.e., periodic changes in the 
fully-polarized ground statepiS, The QD's with different 
impurity locations behave similarly, indicating that the 
circular properties are preserved also when the impurity 
is off-center. This verifies the screening of the impurity 
by the electrons. The importance of screening is evident 
also from the work by Halonen et al^ who considered 
impurity-containing QD's with two and three electrons. 
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The same tendency was also found in the magnetization 
of the two-electron square QD's studied by Sheng and 
Zui^ In our forthcoming studies, we will examine the 
high magnetic field regime as a function of the impurity 
parameters in more detail. 

V. SUMMARY 

We have reported a transport measurement of a verti- 
cal quantum dot, giving single-electron energies strongly 
deviating from the Fock-Darwin spectrum. We have 
found a realistic model for the sample dot by including 
a Coulombic impurity in the external parabolic poten- 
tial that has soft edges. The parameters obtained for 
the model potential indicate a strong confinement result- 
ing from the grown-induced formation of the QD, and 
a repulsive, multiple charged impurity particle shifted 
both laterally and vertically from the QD center. Us- 
ing the model and the spin-density-functional theory, we 
have studied the energetics and structural properties of 



the corresponding many-electron problem. The impu- 
rity evens out the state alternation as a function of the 
magnetic field by lifting the degeneracies, but it does not 
affect remarkably the onset of the full spin polarization. 
The stability of both the maximum-density droplet and 
the magnetization verify that the many-electron struc- 
ture reduces the distortion apparent in the single-electron 
properties. The future experiments may illustrate the 
damping of the irregular behavior in actual quantum-dot 
devices as the number of electrons increases. 
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